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1. Abstract 


The Madhava-Gregory series which is used in the Leibniz formula for z is said to converge 
slowly. This notion deems the formula less useful for practical applications rendering other 
formulas more prominent. But this same formula with no modification but a little change in 
perspective can make it quite competent with others. This paper shows how to improve the 
convergence of the series and make the formula fit for practical applications. Additionally, 
this paper also shows that a basic trigonometric limit identity, combined with an infinite 
precision in calculations, using basic arithmetic can be used to compute the value of z. 


2. Introduction 


During an attempt to build an experimental mathematical system, | needed to use an 
existing mathematical system that provided proofs for the relationships between the 
various formulas in the experimental system. | chose geometry and basic trigonometry to 
address the proofs that were needed. But | encountered an anomaly in the formulas of the 
experimental system when | tried to use them to derive the value of 7. However, during this 
time, | learnt that the Madhava-Gregory series could be easily used to get the value of z in 
other ways than just using a single case where x = 1, which is nothing but the Leibniz 
formula for calculating value of 4. While searching the internet, | was not able to find 
enough information on the matter of increasing the convergence of this series. 


| believe that | have to justify the need to publish such a paper. And so, | wish to mention 
that while experimenting with various methods in the experimental system, | noticed that in 
whatever approaches | used, all values would cancel out, and eventually the Madhava- 
Gregory series would always turn up as the final equation for a. This resulted in instances 
where | could see a pattern in the change of convergence. But this information was not 
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available in any documents that | searched on the internet. Besides, at this point (at the time 
of writing of this document), all information that | can collect on the Madhava-Gregory 
series points to its inability to converge fast enough. So, | have proceeded to derive a proof, 
even though it is quite superfluous, and explain how to work with the convergence of this 
series. It should be noted that, though this proof is an obvious relationship, | did not arrive 
at the proof from the method | show in this document. Actually, | found solutions to various 
cases which resulted in the Madhava-Gregory series, and | saw a pattern from which | did an 
inductive proof. But the one | give here is a simple and straightforward mathematical 
relationship. 


During the development of the formula to increase the convergence | also came up with a 
way to practically show realization of the trigonometric limit identity of x tan(<). Using this, 
| show here that we can get a rudimentary algorithm that uses simple arithmetic to calculate 
the value of z. 


3. The Madhava-Gregory series for 7 


The Leibniz Formula for z, which is actually widely known in the form of the value for 7, is 
an infinite series for arctan(t) for ¢ = 1 and is given as follows: 


114141 
die vere le 1 
“eyo a (1) 


The underlying series is nothing but the infinite Taylor series of the inverse tangent function 
with the expansion at the origin, and is given by: 


arctan(t) =t— —+—>- >+—-—4+... (2) 


which, in sigma notation is given as: 


arctan(t) = 3 (F1)" anes (3) 


a 2n-+1 


When a value of 1 is plugged into the variable ¢ the function yields a value of 7. 


13 1° 17 19 yi 
t 1)=1-—- —+ — —- — + — —- — 4... 4 
eat) a. s eo Ss (4) 


In trigonometry it is known by the integral: 


https://eplogica.com/teplics/mgl-series/icmgl-01/ 2/13 


3/15/23, 4:00 PM Improving the Convergence of Madhava-Gregory Series and a Rudimentary Calculation for tr 


arctan(t) = | : dx (5) 


1+ 2? 


This formula is a combination of Cavalieri's quadrature formula with Maclaurin series for 


arctan(z) !"! that gives the area under the arc of a circle and the series converges for || < 1. 
One specific statement from the internet says: 


Because of its exceedingly slow convergence (it takes five billion terms to 
obtain 10 correct decimal digits), the Leibniz formula is not a very effective 
practical method for computing iT. Finding ways to get around this slow 


convergence has been a subject of great mathematical interest. |! 


This issue of convergence is already addressed in the series itself. The idea is to see what 
affects the convergence. It can be noticed that to increase convergence the value of t needs 
to decrease. Then all we need to do is find values of t < 1 that have a connection with z. 


We already know what happens when t = 1 in (2). And we can easily verify the following 
using the series: 


1 
‘z= arctan(—z). 


e 4 =arctan(2— 73). 


e And go on... 


We can also see it diverging when you have: 
e = =arctan(V3) 


Therefore, we can conclude that for a result = when z is 4 or greater, we can get the series 
to converge. Now you begin to see a way to prove it by induction. However, a simpler proof 
exists. 


4. The proof for the value of z from the Madhava- 
Gregory series 


We know the following from the definition of tan that if: 


tand=t (6) 


then: 
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6 = arctan(t) (7) 
From the understanding of inverse of trigonometric functions we can say: 


6 = arctan(tan 0) (8) 
Here, we can make @ a fraction of 7 by saying: 


= = arctan(tan =) (9) 


We now bring in the Madhava-Gregory series (we will use the sigma notation): 


Ss (=1)" ont 
arctan(t) = as (10) 
dy antl 


and replace arctan(t) and ¢ with values from (9). What we get is: 


T us = {=))" W \on+1 
z= arctan(tan re = ) (tan ra (11) 
Rearranging we get: 


(where K > 4) (12) 


The equation (12) could become the most elegant computation for the value of 7 and the 
proof could be complete. But there is one issue. There is a 7 inside the R.H.S of the 
equation. We may not be able to write a clean equation without 7 on both sides - but we 
can look at K for the solution. However, before we address the issue, we will note a few 
points. They are as follows: 


e In the equation (12) the value of K becomes our constant. 


e Increase in K means decrease in the value of tan — which in turn increases the 
convergence of the series. 


e [tis important to note that we can't use just any arbitrary small value in place of tan <. 
Because it must be linked with K. This is because we need K to be ‘multiplied with' 
after the series is summed up. 


e We can say T = tan =, making T another interlinked constant. 


e Both K and T need to be obtained and they need to be linked. 
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e In practical applications obtaining K is a challenge. This is because from what | could 
collect, there may not be a direct formula to calculate the value of tan —. What we 


have is one for tant!"! which is a series in itself: 


t®? 024) 17t7 ~=— 622? 
tan(t) =t+ —+—-+ + 


aaa (13) 
3 15 315-2835 


But then we have to use the value of z to begin with, and insert it into the formula to 
get the value of z. So, we are back to square one. 


e We don't need a specific value of t, we only need a value of t such that t is very tiny, 
and the associated angle has a form +, and we need both K and T. 


5. Algorithm for z - Obtaining K and T 


There may be a more efficient method to get the value of K by using a straightforward 
mathematical equation (but | am unaware of it). | have tried various methods, unfortunately, 
all ending up with nested square roots. So, we need to deal with nested square roots. But 
because square roots take up a lot of processing time, we need to minimize the operations 
to get an efficient algorithm. In this paper | will provide an iterative programmatic method | 
used, during the value of x proofing experiments, to obtain K with basic mathematical 
operations. 


The idea is to select one angle that is a ratio of a for which we know the value of its inverse 
trigonometric function. With this as our starting point we keep dividing the angle using half 
angle formulas. If we keep track of the number of times we divide the angle, that will 
contribute to the K, in the course we get the value of T as well. Then we use these as 
constants in (12) and choose the number of series iterations to suit our precision of 7. Here 
is how we do it: 


We know that the two trigonometric functions of a fundamental commonly used angle = 
results in values as follows: 


(14) 
(15) 


We are also aware of the trigonometric half-angle identities for cos 8 and tan @ as follows: 
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0 1 — cos 

ail’ 7 / 1+ cos 0 (16) 
0 1+ cosé 

cos = oe (17) 


We see that cos connects with tan very well in the half angle suite. We can now use 
cos($) = + as a Starting point and create an algorithm and iteratively divide the angle as 


we go. Our equation from (12) looks as follows: 


(18) 


We need the values of kK and T so we can plug it into the series to get P which is the value 
of x. R +1 will be the number of times we divide the angle. This we term as the TAN-DIP. L 
will be then number of iterations we run on the series. The algorithm goes as follows: 
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Start TAN-DIP iteration with initial value of 1/2 


ine (oe (19) 


After the TAN-DIP iterations, get the value of tan 


jogo (20) 
1l+vpR 


Determine K 


K=3x2*" (21) 


Put the value into the series 
A (ny 
Pak >. (22) 


Following are some of the observations on this algorithm: 


os 
v2 
K =4 x 2®*1, But we need to make sure that we do at least one of the cos half-angle 


e We can choose the starting point from = cos 7. Then we have to change (21) as 


iterations before we do the tan half-angle calculation. Any starting point below 
1 = cos = can be used with the right initial numbers. 


e We cannot choose starting point from 1 = cos +. Which is obvious. 
e If L = oo and we have infinite precision the P = za. 


e K is proportional to 2”*! therefore increasing the TAN-DIP increases the convergence. 
In other words, higher the TAN-DIP, smaller the angle, greater the K, smaller the T, 
steeper the convergence. 
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e Once K and T are calculated for a very high precision and a large TAN-DIP, it takes 
very few iterations for the Madhava-Gregory series to converge to admirable precision. 


e TAN-DIP iterations are more processor intensive because of the square root 
calculations in each iteration. 


e The algorithm is in its pure form without any numerical adjustments, directly connected 
to the basic trigonometric identities and power-series functions. 


e Any adjustments visible in this algorithm like precision and iteration control are for the 
sake of practicalities. 


6. The Rudimentary Algorithm for z 


There is one more special observation that is worthy of a separate mention. If we take the 
algorithm at (22) and increase the TAN-DIP to a very high value, the convergence goes up 
and it helps us to use lesser and lesser iterations of the series to get the value of 7. We can 
even reduce the iteration to just 1 and still get a decent number of digits to match. This 


stems from the trigonometric limit identity!>! given as: 


lim (n tan =) =T (23) 
n 


n—-oco 


This can be seen not just on paper but in programming reality and it is amazing to witness. 
In the above algorithm: 


e K in (21) is the n of (23) 
e T in (20) is the tan = of (23) 
Then equation (23) becomes: 


lim (KT)=7 (24) 


K-0o 


So if we get very super-high precision, super-high TAN-DIP K and T and just multiply as 
below we get a decent value of z: 


P=KT (25) 


This makes the whole algorithm shrink to: 
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a 1+vR-1 
a 2 
_ l—vpR 
= 1+ up 
K=3~x 2+} 
P= KT 


1! NO SERIES REQUIRED !!! 
Following are some of the observations on this algorithm: 
e K is proportional to 2#+! therefore increasing the TAN-DIP increases the precision of 7 
digits. 
e The algorithm is in its purely arithmetic form and without any numeric adjustments. 


e Related formulas which are close to this one are Viéte's formula!® and Liu Hui's 7 


algorithm !71, But both of them are of geometric or geometric-trigonometry origin. This 
algorithm is of the trigonometric-calculus origin. 


7. Experiments and Results 


A simple code with low TAN-DIP, is written in PHP and shown below. With low TAN-DIP it 
doesn't amount to much. TAN-DIP calculations take a lot of time as well. 


g$iterations = 5; 
$pi_decimals = 2000; 
$precision = $pi_decimals*3; 


bcscale($precision) ; 
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// TAN DIP 

$v = bcdiv(1,2); 

$j = $R - 1; 

for ($i=0; $i < $j; $i++) { 
$v = bcsqrt( 

bcdiv(bcadd(1, $v), 2) 

)3 

} 


$K = bcmul(3, bcpow(2, $R)); 
$T = bcsqrt(bcdiv(bcsub(1, $v), bcadd(1,$v))); 


// THE SERIES 
g$sum = @; 
for ($i =0; $i < $iterations; $i++) { 
g$exp = 1 + (2 * $i); 
$sum = bcadd( 
$sum, 
bcdiv( 
bcmul(pow(-1, $i), bcpow($T, $exp)), 
$exp 


)5 
} 
$P = bcmul($K, $sum) ; 


Snapshots of the results are shown below, they are very modest because the convergence is 
slow due to low TAN-DIP and the total series iterations are controlled to 5. It managed to 
match 31 digits. 


actual 3.141592653589793238462643 383279502884 1971693993751 05820974 
matched 3.141592653589793238462643383279860 194 12395 7232143013882393 
matched_digits 31 

precision 6000 

iterations 5 

R 10 

T 0.001022654215094841154940098025921773462183801600676542003 
K 3072 

tandip-time 15978mS 

series-time 452mS 


Figure 1: Simple Algo Sample Result 
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During the development for this paper multiple experiments were carried out with various 
precisions and TAN-DIPs. These were done as a test to see what the effects were of various 
parameters like K, T, L, R etc. But they were not designed to profile or measure the 
performance as it was not a part of the study. During these experiments the following was 
observed: 


e TAN-DIP needs to be at least 1000 to get good admirable results and the calculations 
must be done with at least 3 times the precision required. 


e On many experiments TAN-DIP of 1001 with precision 6000 (K-T-6000-1001) was used 
as a constant. 


e With K-T-6000-1001 as a constant, 7-precision even up to 3000 digits could be got 
with 5 iterations. 


e Some examples are shown below: 


actual 3.1415926535897932384626433832 79502884197 1693993 75105820974 
matched 3.141592653589793238462643383279502884197169399375105820974 
matched _digits 3014 

iterations 5 

R 1001 

Vv 0.00000000000000000000000000000000000000000000000000000000C 
K 642905 1643 117603925690550294360010863368428870233201644662: 
series-time 449mS 


Figure 2: With large R and high precision but pre calculated constant. 


actual 3,141592693589793238462643383279502884197169399375105820974 
matched 3.141592653589793238462643383279502884197169399375105820974 
matched_digits 603 


iterations 0 

R 1001 

Vv 0.00000000000000000000000000000000000000000000000000000000¢ 
K 6429051643 1176039256905502943600108633684288 70233201644662: 
algo 2mS 


Figure 3: An example of just multiplying K with T with large R and high precision. 


actual 3.1415926535897932384626433832 79502884197 1693993 75105820974 
matched 3.141592653589793238462643383279502884197169399375105820974 
matched_digits 2037 

iterations 146 

algo 2969mS 


Figure 4: An example of Chudnovsky's algorithm based results 
(Algorithm - Courtesy of https://github.com/natmchugh/pi) 
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Some points to note: 


e This paper deals with the proofs of the Madhava-Gregory series and how to deal with 
convergence and get results but not the quality of results. 


e The various practical aspects such as profiling, and performance of this algorithm are 
not the focus of the paper. The experiments to study the effects of variations in the 
parameters such as K, T, L, R etc. must be done separately. 


e All experiments were done in PHP 7.3.4 on a Windows 1@ PC with Intel(R) Core(TM) 
i7-4770 CPU @ 3.4@GHz and 24.@ GB RAM. 


8. Conclusion 


These conclusions are voiced as opinions and no more: 


e The Madhava-Gregory series is a natural (free of numeric adjustments) method to 
calculate the value of z. 


e Moreover, because this formula has its basis in straightforward simple trigonometry, it 
can be the choice to provide the simplest of proofs for the value of z. 


e In terms of the algorithm we have seen that convergence could be improved using (K, 
T) which is an interlinked constant pair. 


e The profiling, performance, radius of convergence, parametric analysis etc., of the 
algorithm must be a part of another study. 


e If calculation of (K,T) is not considered as part of this algorithm, then this algorithm, 
with a good (K,T), performs superior to the existing Chudnovsky’s algorithm. But such 
a claim would be wrong as the pair (K,7) itself is part of the calculation and the time 
taken by the algorithm, because the pair is part of the calculation of 7. Chudnovsky's 
algorithm is currently still the best at speed. 


e In terms of the Spirit of Mathematics, just as we have rudimentary geometric 
calculations for the value of a through iterations of polygons, we also see that a 
rudimentary calculation for z in terms of pure arithmetic exists through the knowledge 
of the trigonometric-calculus. 
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